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We investigate the von Neumann entanglement entropy as function of the size of a subsystem 
for permutation invariant ground states in models with finite number of states per site, e.g., in 
quantum spin models. We demonstrate that the entanglement entropy of n sites in a system of 
length L generically grows as a log 2 [2nen(L — n)/L] + C, where a is the on-site spin and C is a 
function depending only on magnetization. 
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I. INTRODUCTION 



Recently it has been argued that for critical (gap- 
less) quantum spin systems the entanglement entropy of 
a block of n spins diverges logarithmically as 7log 2 ?i, 
while for non-critical systems it converges to a constant 
finite value 0, H 13. This property was interpreted in 
the framework of conformal field theory Q associated 
with the corresponding quantum phase transition and 
the prefactor 7 of the logarithm was related to the central 
charge c of the theory as c = 37. This was shown explic- 
itly for the exactly solvable antiferromagnetic Heisenberg 
spin 1/2 chain, i.e., the XX Z model where the anisotropy 
parameter A belongs in the critical regime to the interval 
(—1,1) and 7 = 1/3. Rather surprisingly, at the transi- 
tion point A = — 1 from gapless to noncritical behaviour, 
the entanglement (von Neumann) entropy of a block of 
spins in the ground state was found to grow faster than 
in the critical domain — 1 < A < 1, namely as 7log 9 n 
with the logarithmic prefactor satisfying | < 7 < 1 [5j. 
At A = — 1 the ground state of the XX Z Hamiltonian 
has permutational invariance (up to a gauge transfor- 
mation), and is degenerate with respect to the total z- 
magnetization, so that the whole system is generically 
described via a density matrix which can be written as 
a (weighted) sum of projectors on the multiplet compo- 
nents. The lower bound 7 = i is attained if the state of 
the whole system is a pure state (only one component of 
the multiplet is present), while the upper bound 7=1 
is reached for a mixed state where all the components 
of the multiplet have equal weights. Note that the pure 
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state is generically not factorizable, see (@J below, thus 
producing a mixed state after partial tracing. 

In the present paper we generalize the approach of || 
to arbitrary permutation invariant quantum spin states. 
In particular, we consider the case of a ferromagnetic spin 
chain with arbitrary spin a on every site and show that 
the entanglement entropy of a block of n spins in the fer- 
romagnetic ground state generically diverges as a log 2 n. 
Our approach uses the invariance of the ground state un- 
der site permutations, thus allowing us to compute the 
entanglement entropy exactly for blocks of arbitrary size 
and for systems of arbitrary length. 

The paper is organized as follows. In section 2 we in- 
troduce the permutation invariant states and list physical 
systems whose ground states have this symmetry. In sec- 
tion 3 we formulate a theorem which gives the analytical 
expression of the eigenvalues of the reduced density ma- 
trix for arbitrary spin a. Using this theorem we compute 
the entanglement entropy of a block of size n in the finite 
system of total length L in specific ground state sectors. 
Taking the limit of large subsystem sizes, we derive ana- 
lytic expressions for the entanglement entropy Sr n ) both 
for n,L»l and for n^$> 1, L = 00. As a result, we find 
that in the ground state sector with a fixed value of S z 
the block entanglement entropy diverges for large n as 
Sf n ) = fjlog 2 [27ren(L — n)/L] + C. In section 4 the case 
of spin cr > 1 is treated in more detail. A discussion and 
some further remarks close the paper. 



II. PERMUTATION INVARIANT STATES 

Let us consider states in a Hilbert space (C d )® L of a 
quantum spin chain of local spin a, where d = 2a + 1 
is the dimension of the spin space at every site i and d L 
the dimension of the whole Hilbert space. Permutation 
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invariant states constitute a subspace Q of substantially 
smaller dimension 



k{L) = dimQ = 



L + 2a 
2a 



(1) 



this being the number of possible ways to distribute L 
indistinguishable objects among 2a + 1 boxes. Denoting 
by Nj the number of objects in the block j, the state of 
the whole system is completely characterized by 



\9(L,N ,N u ...,Naa)) 



N \Ni\...N 2 J 



(2) 



N, 



N 2 „ 



where iVo is the number of spins pointing down (a z = 
— cr), Nj is the number of spins with a z = —a + j , up 
to N 2a spins with maximal a z — a, occupying in total 
L = No + Ni + ... + N2a sites. The sum is taken over 
all possible permutations, the total number of which is 
n \Ni\ n-2 ! ' w ith the prefactor in (J2J taking care of the 
normalization. 

We will be interested in the entanglement of a block 
of n spins with the remaining L — n spins (playing here 
the role of the environment), characterized by the von 
Neumann entropy 



S( n ) = -tr(p n log 2 p n ) 



k 



X k log 2 Afe, 



(3) 



where p n is the reduced density matrix of the block, 
obtained from the density matrix p of the whole sys- 
tem by tracing out the degrees of freedom of the en- 
vironment p/ n \ = £r (£_„)/?, and Xk are its eigen- 
values. The density matrix of the whole system is 
a projector on the pure state in (J2J), i.e. p = 
\y{L,N ,N 1 ,...N 2 < T ))(*(L,N ,N 1 ,...N 2<T )\. Notice that 
due to the permutational symmetry, SV n ) does not de- 
pend on the particular choice of the sites in the block 
but only on its size n. The eigenvalues Afe of the reduced 
density matrix are by construction all real, nonnegative, 
and sum up to one: Xk = 1. 

Before giving the general expression for Afe we discuss 
separately the two-states case of a = 1/2. In this case 
the state of the system |(2J reduces to 



\*(L,N)) 



-1/2 



X) in -ni- 1) (4) 



TV 



L-N 



where f and J, denote spin up and spin down respectively. 
This state appears in the literature in different physical 
situations. Since the entanglement properties do not de- 
pend on the underlying model but only on the form of 
the wavefunction, we list some models for which Q is an 
exact ground state to show the diversity of applications. 



i) The isotropic Heisenberg ferromagnet, H = 

J Yli=i Vi^i+i where ai are Pauli matrices, J < de- 
notes the exchange constant and L is the total num- 
ber of spins. (We assume periodic boundary conditions 
L + 1 = 1.) The ground state belongs to a multiplet of 
total spin S = L/2 and is (L + l)-fold degenerate with 
S z = + 1, The state Q is a pure state, 

corresponding to the multiplet component with a fixed 
number N of spins up. 

ii) The Heisenberg antiferromagnet at A = — 1. In this 
case, the state (@J is obtained from the ground state by a 
unitary transformation, inverting every other spin along 
the chain (note that the von Neumann entropy is invari- 
ant under unitary transformations of the state of the sys- 
tem). 

iii) The generalized Hubbard model in the limit of strong 
attraction (the so-called eta-pairing stateslg). 

iv) Hardcore bosons on a complete graph [J, 

v) The Lipkin-Meshkov-Glick model [Tl|. 

The state Q also appears in a classical context as sta- 
tionary distribution of the asymmetric exclusion process 
(ASEP) describing a Markov process of noncquilibrium 
stochastic motion of N particles on a ring with L sites 
with hardcore exclusion. In this case (@J means that all 
particle configurations have equal probabilities in the sta- 
tionary state 01 • The block entanglement entropy for 
the case of spin cr = 1/2 was obtained in 0. In particu- 
lar, it was shown that the eigenvalues Afe of the reduced 
density matrix p n are Xk{L, n) — CjtC^ZJk /Cjv > where 
denotes the binomial coefficient n\j (m\(n — m)!) and 
k = 0, 1, ... min(n, N). In the limit of large n, the von- 
Neumann entropy was found to be 

1 1 n(L — n) 
S(n)(P) ~ g lo S2W) + -log 2 27re . (5) 

where p=^,q=l—p. In the next section we generalize 
these results to the case of arbitrary spins. 



III. ENTANGLEMENT ENTROPY FOR 
ARBITRARY ON-SITE SPIN a 

We formulate the main result of this section in the 
following 

Theorem: The eigenvalues of the reduced density matrix 
p n (No, Ni, ...N2a) of a block of n spins in the permuta- 
tion invariant state 0) of the whole system are given by 



2a 

X k (L,n,a) = - I l[c£ 



(6) 



where ki,N satisfy ko + ... + k 2rT = n, N + ... + N 2rJ = L. 
Here and below we denote the set of ki by the bold k. 
To prove the theorem we note that the reduced density 
matrix p n is decomposed with respect to symmetric or- 
thogonal subspaces of the system of n spins, classified 
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FIG. 1: Exact (points) versus analytical (curves) comparison 
of von-Neumann entropy S( n ) for spin a = 1 in systems with 
L = 30, 60, 120, 240, oo. Blocks of all sizes < n < L (for 
finite L) are considered. The other parameters = = g . 



by the numbers kj = 0, 1, ...min(n, Nj) of spins with 
u z = —a + j in the block 



l (JV)=5^to k |^(n 1 k))(^(n,k)|. 



(7) 



Here the state \ip{n, k)) has the same structure as 



|V(n,k)) = £ | U ..4- ... /" -TT - T) (8) 

and w/c is the corresponding probability, given by the 
number of ways one can distribute Nj spins (of values 
—cr + j) on kj sites for all possible values of Nj,kj, 
j = 0, ...,2a, divided by the number of ways one can 
distribute N = No + ... + N 2a spins on n = fc + ... + k 2a 
sites (the total number of states), i.e. 



1 2<T 



(9) 



Since the states |^>(n, k)) are orthogonal, the represen- 
tation in Q is diagonal and the coincide with the 
eigenvalues of p n (N). This proves of the theorem. No- 
tice that for the case a — 1/2 Eq. © reproduces the 
results derived in p| (use N = L — N\ , fc = n — k\ and 
the invariance under the exchange of N% and n) . Having 
found the eigenvalues of p n (N) one can compute the en- 
tanglement entropy St n ) f° r arbitrary L, n and N . In the 
large-n limit one obtains, see the Appendix, in analogy 



FIG. 2: Entanglement entropy as function of a number of 
sites n in the block, for different on-site spin a = |, 1, |,2, 
for finite system L — 120. Comparison of the exact formula 
(points) with the limiting expression IjlOfl (continuous curves) . 
In all cases equal partitions of particles is taken pi = Ni/L = 

1/(2(7+1). 



with the case of a = 1/2 5], that 
S (n) (L,a) = C + a log 2 (tore "^ ~ n) 

C = \\og 2 (jjpj , 

Vfe=0 / 



n > 1, 
(10) 

(11) 



where pk = In Figs. H 121 we compare the exact en- 
tropy of finite systems, as computed from the exact ex- 
pressions Eqs. with the analytical expression l|10fl . 
from which we see that there is an excellent agreement 



2(T 



also for small values of nJ^Pfe. In the thermodynamic 



k=0 



limit L — > oo, -^j- — > p,-, the eigenvalues of the reduced 
density matrix JoJ) simplify to 



A k (co,n, { Pl }to)^ n\]J 



i=0 



(Pi)" 

hi 



The entanglement entropy then becomes 

S ( n) (oo, fe}?=o) = C{{ Pi }) + a log 2 (2^en) (12) 

where the constant C is given by (|11|) . 

Notice that while the general formulae [lilGII are valid 
for an arbitrary choice of parameters, the analytical re- 

2(7 

suit (|10|1 is valid for nj^pfc ^> 1. If in the limit of an 

k=0 

infinite system L — > oo some of the pk vanish (say, z 
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coefficients pk — > out of a total number of coefficients 
2a), the entangled state of n spins will effectively be- 
have like the one with effective site spin a = a — #, 
and the von Neumann entropy will respectively grow as 
S n ~ (c — |) log 2 ft instead of 5 n ~ u log 2 n, see (|T2*)) . 



It is instructive to recall that Eq. (|10|l is derived un- 
der the assumption of a pure state of the global system. 
If the state of the whole system is mixed in the ensem- 
ble of states with the same symmetry the resulting 
von-Neumann entropy S( n ) will become larger. An upper 
bound of S7 n ) can be derived by noting that the number 
of nonzero eigenvalues of the reduced density matrix is 
equal to the number of terms in the decomposition Q, 
which is «(n), see Q. The upper limit of the entangle- 
ment entropy is thus reached when all in © are equal, 
implying 



sup(5( T . 



log 2 



n + 2cr 
2a 



(13) 



This limit is attained for a mixed state of the global 
system with all components of the ensemble © equally 
weighted. In this case the von-Neumann entropy is given 
by the above expression even for finite n, L (see |5j for an 
example). For large n, sup(SV n \) w 2a\og 2 n. 

We also remark that in the critical spin models where 
spin-spin correlation decay algebraically in the ground 
state (e.g., the region — 1 < A < 1 of the antiferromag- 
netic Heisenberg chain) , there are three distinct physical 
properties contributing to the entanglement: (a) on-site 
correlations due to single-site quantum fluctuations, (b) 
algebraically decaying spin-spin correlations which sur- 
vive in the thermodynamic limit and (c) the correlations 
due to the constraint of fixed magnetization which van- 
ish in every domain of finite size in the thermodynamic 
limit. For permutation invariant states @ only contri- 
bution (a) is left in the thermodynamic limit for finite 
domains, but for domains of finite volume fraction also 
the correlations (c) due to the constraint of fixed magne- 
tization remain relevant. 



IV. MODELS WITH HIGHER SPIN 

As remarked above the entanglement entropy of a 
quantum state does not depend on any underlying model 
but only on the properties of that state. Nevertheless it 
is of interest to have some insight for which systems the 
permutation invariant states J2J considered here are the 
ground state of that quantum system. First of all, the 
generalizations of the Sutherland model ^3| describing 
quantum spin chains with an interaction given by the 
permutation operator in SU(N), 



H 



"A?-^ 



example being nearest neighbor interaction. For ferro- 
magnetic interaction, all < 0, the states span its 
ground state, which can be proved along the lines of [l4| . 
For SU(2) and Jjj = J<5; J+ i the Hamiltonian l|14f> re- 
duces to the isotropic Heisenberg Hamiltonian. For the 
general case we recall that the permutation operator is 
written in terms of spin operators as 



2S 2S2 SOS )-k(k + l) + 2S(S + 1) 

p=£(-i)**n- 

8=0 k^i 



i(i + 1) - k(k + 1) 



where N = 2S + 1, (see, e.g., ITJj)- For the above exam- 
ple, all states © are eigenstates of (|14|l with the lowest 
energy. 

The increase of entanglement entropy compared to the 
s = 1/2 case appears to result from the larger local state 
space of the SU(N) chain, but not from the SU(N)- 
symmetry itself. This picture is supported by the gen- 
eralized disordered SU(2)-symmetric spin-cr Heisenberg 
ferromagnet 



H 



^2 Ji >o9{viVj 

i,3 



1) 



(15) 



where the exchange energies Jjj < may be non-zero 
between any pair of lattice sites on an arbitrary lattice 
and <7i are local SU(2) generators in the spin-cr repre- 
sentation. For any polynomial function g with positive 
expansion coefficients linear combinations of the permu- 
tation invariant states @ with fixed total magnetization 
are ground states of the Hamiltonian, see [16j for a de- 
tailed discussion of the ground states of this model in a 
probabilistic setting. Since these are not pure states in 
the sense discussed above (with all quantum numbers Nk 
fixed) the entanglement entropy of these ground states is 
higher than those of the pure ground states of the SU(N) 
spin chain. 

Other models with pair interaction, but no symmetry: 
Using the Perron- Frobenius theorem it is straightforward 
to construct quantum Hamiltonians of the structure 



H 



(16) 



1,3 



(14) 



where gi d is a hermitian pair interaction matrix satis- 
fying g%j\ s ) = for all i,j and where the wave func- 
tion | s ) which is constant for all spin configurations in 
(<D rf ) 0L is the ground state of g itj . Then, if all J id < 0, 
the vector | s } is also the ground state of if, and if fur- 
thermore gij has no invariant subspaces, it is the unique 
ground state with maximal entanglement entropy l|13fl . 
Such quantum systems also have a probabilistic inter- 
pretation as generator of some irreducible Markov chain 



SUMMARY AND CONCLUSIONS 



are obviously invariant under the permutation group Sn- 
Here the set of J^i is defined on any connected graph, an 



We have obtained exact eigenvalues of the reduced den- 
sity matrix for permutation invariant states, for arbitrary 
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length of the system. In the thermodynamic limit, it was 
shown that the von Neumann entropy of entanglement 
of a block of n spins a with the environment grows log- 
arithmically fast Si n -\ ~ 7 log 2 n , with a prefactor 7 = a 
for a pure global state and with 7 = 2a for homoge- 
nously (maximally) mixed global state. Various models, 
the ground states of which are permutation invariant, 
are given, for spin a — 1/2 and higher. We note that 
the logarithmic growth of entanglement entropy due to 
permutation invariance is faster than the one of criti- 
cal (conformally invariant) models, where S< n \ ~ 7 log 2 n 
with 7 = 1/3 was observed 0. 

It is also interesting to compare the finite size cor- 
rections of the entropy of the permutation invariant 
states lfT7)|) with those of critical spin chains, ~ 

§log 2 (£sin(^)), obtained in \V2} , see Eq.(3.8) in this 
paper. Expanding in the first nonvanishing order of 
1/L both expressions, we get the finite size corrections 
A-(n) = Sfo(L) - S£,(oo), A^(n) = Sf^(L) - 



(«) 



£?^r(oo) for permutation invariant and critical models 



as 



1 n . /nv 
(J ln2Z + (z) 



AP e » = alog 2 (l- z ) 

A-(n) « -log 2 (l - k-f) « f-) 2 + O(-) 

y ' 3 B2V 3 L 9 In 2 V L / \L 

Thus, besides the difference of the coefficient of the prin- 
cipal logarithmic divergence (7 log 2 n with 7 = 1/3 for 
critical XX Z and 7 = a for permutation invariant states 
resp., a substantial difference in the finite size corrections 
(linear in n and of order 1/L for permutation invariance) 
and quadratic in n and of order 1 /L 2 resp. (for conformal 
invariant critical states) is also observed. 
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VI. APPENDIX 

To prove Eq. I|1(J|I we first give the general scheme and 
then demonstrate the result on a particular example. For 
simplicity we shall discuss the limit when the size of the 
global system tends to infinity, L — > 00. The eigenvalues 
(JHJ) Ak = A(fco, fci, ...fc 2o -) for large values of n and fcj, 
can be approximated by a multi-dimensional Gaussian 
distribution with mean and moments given by 



npi 



(17) 



((fcj - (fci)) 2 ) = npi (1 - pi) 



(18) 



{{ki- (ki))(kj - (kj))) = -npiPj, i^j (19) 

where Pi = i = 0, 1, ...2a, is the average density of 
spins —a + i in the system. To prove Eqs. H17 jl -I|19 [l 
we first observe that, using the Stirling approximation 
ml = m m exp{—rri)\j2-Km, in the limit L,N{ — > 00 the 
eigenvalues © of the reduced density matrix can be ap- 
proximated as 



A k ~ Ck(n) Y[pi* 



(20) 



i=0 



with Ck(n) — 7j!/(fco!fci!...fc2 CT !)- Notice that in this limit 
the distribution of the Ak coincides with the multinomial 
distribution 

/ 2(7 \ n 2(T 

(Eft) =E*NlIft fei =E^ = 1 - 



i=0 



Using this expression, the mean value of fc Q , a — 
0, 1, 2a, can be found as 



2a 



d 



(k a ) = E fcaCk (")IT^ fel = Pa gf(H Pi ^ 



k i=0 

In a similar manner we obtain 
d 2 



np 



i=0 



(k a (fc Q - 1)) =P a j^2 ( E K ) = n ( n_1 )Pa> 

( la \ " 

{k a kp) = PaP/3 a „" a „ ( E Pi ) = ~ l )P a PP' 



'dpi 

d 2 



dpadpp 



,i=0 



from which Eqs. I|18|) . I|f 9|l readily follow. Denoting by 
Xi = ki/n, so that < x% < 1, the eigenvalues Ak are 
approximated as 



let A ( 1^ A 

— exp I — - A y (pa - p.,) [Xj - pj) J 



Vdet A 
72" 



where the shifting of introduced to account for 

the nonzero mean l|17(l . and the coefficients of the sym- 
metric matrix <Zjj = a,ji are fixed using the moments 11811 . 
(|19|l . On the other hand, the computation of the mo- 
ments from the distribution (|21|) gives 



{{Xi-Pi) (Xj - pj)) = 



detA 



(22) 



where My- are the minors of the matrix A. Comparing 
Eq. (|2*2*|) with Eqs. lfT%l ). lfT§)l . we can fix the elements A^ 
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of the matrix. To illustrate the method let us consider 
the case of spin a = 1. Then, 



A k = AO, y) » exp (- i (ax 2 + 26a;j/ + cj/ 2 



where we introduced x = X\ — pi, y = X2 — P2 and D = 
detA = ac — b 2 for brevity of notation. From (|18|l . (|19|) 
we have (x 2 ) = npt(l - p\), (y 2 ) = np 2 (l - P2), and 
(xy) = —npip2- On the other hand, from 1221) we obtain 
(x 2 ) = c/D, (y 2 ) = a/D, and (xy) = -b/D. Computing 
the determinant 

D = ac-b 2 = D 2 {(x 2 )(y 2 )-(xy) 2 ) 

and substituting the moments, we get 

D~ x = n 2 pip 2 p3- 

The von Neumann entropy is then computed as 



S 



dx / dy(X(x,y)log 2 X(x,y)) 



Notice that for finite pi the larger contribution to the 
integral comes from the neighbor of the origin x = y = 
so that we can extend the limits of integration to the 
whole real axis and perform the integral exactly. This 
leads to 



S(n) a log. 



2nire 



D 



log 2 2?me + \ log 2 (P1P2P3) , 



which coincides with the expression l|10|) in the limit L = 
00 for the case of spin a = 1. A more detailed analysis, 
analogous to the one done in 0, restores the finite-size 
dependence on L of the von-Neumann entropy as 



S(n) a log 2 2-7re 



n(L — n) 



+ 5 lo & {P1P2P3) for a = 1 



Working out the same procedure for an arbitrary spin a 
we obtain (|10() . 
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